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Abstract. In this paper we will explicitly work out the complex first-order SUSY transfor- 
mation for the harmonic oscillator in order to obtain both real and complex new exactly- 
solvable potentials. Furthermore, we will show that this systems lead us to exact complex 
solutions of the Painleve IV equation with complex parameters. We present some concrete 
examples of such solutions. 
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1 Introduction 

Supersymmetry was introduced in quantum field theory in order to relate the bosonic and 
fermionic sectors into one unified superfield [37]. This theoretical framework is described by 
some deformations of Lie algebras and their commutation relations. Although there is no expe- 
rimental evidence of supersymmetry in nature yet, its concepts have aroused new ideas in other 
fields of science. 

The supersymmetric quantum mechanics (SUSY QM) is one of these new ideas: while first 
developed by Witten [39], it was soon related to the Darboux transformation and to the facto- 
rization method, which was thought to be completed in 1951 by Infold and Hull [25] but it was 
developed further as a paradigm to obtain new exactly-solvable by Mielnik [30] and others [6, 
7, 17, 31, 35] and deepen into the algebraic structure of quantum systems. Even more, for 
those potentials that are not exactly-solvable we still have a powerful new tool to develop 
approximation methods, giving rise to a whole new branch of study called spectral design [5], 
since one can obtain an exactly-solvable potential with a prescribed spectrum, i.e. with some 
given energy levels. 

In this paper we will use the scheme of SUSY QM to obtain complex partner potentials. 
As far as we know, they were first developed in [4], where complex potentials with real energy 
spectra were obtained using complex transformation functions associated with real factorization 
energies and later in [19] with complex factorization energies. On the other hand, in this work 
we will study complex potentials with complex energy spectra. 

The theoretical framework of SUSY QM have connections with other branches of study, for 
example with non-classical orthogonal polynomials and the special functions [34], which are so- 
lutions of linear differential equations. In this paper, we are interested in some other functions 
that play an analogous role as special functions, but now for non-linear differential equations. 
Specifically, we are talking about the Painleve equations. Some specialists (for example [15, 26]) 
consider that solutions of the Painleve equations will be the future members of special functions 
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in the twenty-first century. It turns out that one of these equations, the Painleve IV equa- 
tion (.Piv) (see [24]), is closely related to potentials obtained by SUSY QM from the harmonic 
oscillator. Note that Pry also appears in other areas of physics, for example in fluid mecha- 
nics [38], non- linear optics [21], and quantum gravity [22]. It also arises in mathematics, e.g., in 
the symmetric reduction of several partial differential equations like Boussinesq equation [14], 
dispersive wave equation [32], non-linear cubic Schrddinger equation [12], and auto-dual Yang- 
Mills field equations [1]. 

The SUSY QM technique has been used before to obtain solutions of Pry; however, they 
are either real or complex but associated with a differential equation having real parame- 
ters [9, 10, 11]. In this work we will expand the solution set to obtain complex solutions, 
associated now with complex parameters of Piv- 

To accomplish this, the structure of the paper is the following: in Section 2 we will work out 
the complex scheme of the first-order SUSY QM. In Section 3 we will describe the second-order 
polynomial Heisenberg algebras and its relationship with Prv- Then, in Section 4 we will use the 
first-order SUSY partners of the harmonic oscillator in order to connect with complex solutions 
to Piv- These solutions shall be studied in Section 5. Finally, we will present our conclusions in 
Section 6. 



2 Complex SUSY QM 

Let us begin with a given Hamiltonian H which has been completely solved, i.e., all of its 
eigenvalues E k and eigenfunctions ip k , k = 0,1, . . . are known: 

H = -d 2 + V(x), (1) 
H^ k {x) = E k i/> k (x), (2) 

where d = d/dx and we are using natural units such that h = m = 1. 

Now, we use the factorization method which is equivalent to SUSY QM but much simpler in 
this case. We propose that H is factorizable as 

H = A'A + + e. (3) 

In the standard real factorization method there is an extra condition A + = (A~)^ . In this paper 
we do not use this constrain, but rather we simply ask that 

A+ = -d + f3(x), A- = d + (3(x), (4) 

where @{x) is a complex function to be found. This choice represents a more general factorization 
than the usual real one [33] . 

Working out the operations in (3), using the definitions in (1) and (4), we obtain one condition 

for j3(x) 

f3' + f3 2 = V(x) - e, 

which is a Riccati equation. 

On the other hand, if we consider a similar factorization but in a reversed order and introduce 
a new Hamiltonian H, defined by H = —d 2 + V and 

H = A + A-+e, (5) 

it turns out that 



V{x) = V{x) - 2f3'{x). 
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Besides, from equations (3) and (5) it is straightforward to show that [4, 33] 

HA + = A + H, HA~ = A~H, (6) 

which are the well known intertwining relationships with A + , A~ being the intertwining opera- 
tors. From equations (2) and (6) we can obtain the eigenvalues and eigenfunctions of the new 
Hamiltonian H as follows 

HA+^ k (x) = A + H^ k (x) = E k A + M*), H [A+M*)} = E k [A + ^ k {x)] . 

Therefore, the eigenfunctions ip k of H associated with the eigenvalues E k become 

ip k oc A + tp k (x) oc W ^ k , 
u 

where W is the Wronskian, (3(x) = [In u(x)]' , and u(x) is an eigenfunction of H (in general 
a non-physical one) associated with a complex eigenvalue e, i.e. 

Hu = eu. 

Furthermore, the eigenstates ip k are not automatically normalized as in the real SUSY QM since 



now 



(A + i> n \A + il, n ) = (if> n \(A + ^A + ip n ), 

and in this case (t4 + )^A + ^ (H — e). Nevertheless, since they are normalizable we can introduce 
a normalizing constant C n , chosen for simplicity as C n G M + , so that 

i> n (x) = C n A + Tp n (x), (i> n \lpn) = 1- 
Finally, there is a wavefunction 

u 

that is also eigenfunction of H 
Hip € = eip e . 

If it is normalizable, it turns out that V(x) is a complex potential which Hamiltonian H has the 
following spectrum 

Sp(tf) = {e}U{£ n ,n = 0,l,...}, (7) 
and e G C although in particular e could be real. 

3 Polynomial Heisenberg algebras and Piy 

A polynomial Heisenberg algebra (PHA) is a particular deformation of the Heisenberg-Weyl 
algebra [18, 36]. Although there is a general definition of a m-th order PHA, in this work we 
will only be interested in the second-order case. 
A second-order PHA is defined by 

[H,^] =±2L ± > [L~, L + ] = Q(H + 2) - Q(H) = P(H), 

where P and Q are second- and third-order polynomials, respectively. 
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Figure 1. Diagram of the two equivalent SUSY transformations. Above: the three-step first-order 
operators Lf, L 2 , and L 3 allow to accomplish the transformation. Below: the direct transformation 
achieved through the third-order operators L , 

This definition means that act like ladder operators for the system described by the 
Hamiltonian H [13, 20]. As usual, L + will be a creation operator and L~ an annihilation one, i.e., 

In this case, both will be third-order differential ladder operators. Now, we propose 
a closed-chain of three SUSY transformations [2, 3, 16, 23, 28, 36] so that can be expressed as 

L+ = L 3 +L+L+ = (d - h)(d - f 2 )(d - h), 

L~ = L^L 2 L 3 = (-d - h){-d - f 2 )(-d - h). (8) 

In general, {L~)^ 7^ L + , except in the case where all fi GR. As a matter of fact, in this article we 
are interested precisely in the case where fi R. Each Lj fulfills two intertwining relationships 
of kind 

Hi+iLf = LfHi, HiLj = L i fl»+i, (9) 

where i = 1, 2,3. In Fig. 1 we present a diagram of the SUSY transformations. 

If we equate the two different factorizations associated with (9) which lead to the same 
Hamiltonians we get 

Hi+i = + Si = L i+l L+ +l + ei+i, i = 1,2. 

In addition, the closure condition is given by 

H4 = + €3 = H% — 2 = + 61 — 2. 

By making the corresponding operator products we get the following system [2, 29, 36] 

(10a) 
(10b) 
(10c) 

Eliminating /| from (10a) and (10b) we get 
fi + V'i + f'^fl -fl + ei-e 3 , 
and from here we substitute f\ from (10c) to obtain 

/l + /2 + /3 = l> 

which, after integration becomes 

h + h + h = x. (11) 



/1+/2 


= /i 2 


-/! 


+ ei 


- £2, 


fL + fa 


= /! 


-/I 


+ e 2 


- £3, 


/1 + /3 


= /I 


-A 2 


+ e 3 


-ei + 2 
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Now, substituting (11) into (10a) 

x - /3 , I-/3 ei - e 2 



2 2{x-h) 2(x-f 3 Y 
Let us define now a useful new function as g = fa — x, from which we get 
9 9' ei -62 

Similarly, by plugging (11) into (10b) and using 5 we obtain 

, _ _9 _ (/_ _ gi ~ e 2 
2 2 5 2 5 ' 

Now that we have /1, /2, fs in terms of 5, we replace them in (10c) in order to obtain 

" 1 I '\2 , 3 4 , a 3 , n 2 f 2 , , i 61+62^ (ei — 6 2 )" 

55 =2^) + 2 9 + 4 5' a; + 2 5 f (a: + d + 1 ^- ' 

which is the Painleve IV equation (Piy) 

99" = \{g'? + + + V (x 2 -a)+b, 

with parameters 

ei + e 2 (ei - e 2 ) 2 
a = -2--e 8 -l 1 6 = . 

Since, in general / E C then g £ C. In addition, £j S C which implies that a,i 6 C and so g is 
a complex solution to Prv associated with the complex parameters a, b. 



4 First-order SUSY partners of the harmonic oscillator 

In this section we will show that the complex first-order SUSY transformation applied to the 
harmonic oscillator leads to a system with third-order ladder operators which is naturally ruled 
by a second-order PHA. 

Let us begin with the harmonic oscillator potential 

V{x) =x 2 . 

As in Section 2, we propose two operators A + and A~ that fulfill equations (4) and (6) so that 

V = V-2(3'(x), 
with the condition that j3{x) should solve the Riccati equation 

p + /3 2 = V - e. (12) 

The previous scheme can be written in terms of a solution of the initial stationary Schrodinger 
equation u by substituting f3 = (In u)' in (12) in order to obtain 

— u" + Vu = eu, 
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Figure 2. Examples of SUSY partner potentials of the harmonic oscillator using the two complex 
factorization energies e — + i with A = re = 1 (left) and e = 3 + ilO -3 with A — re = 2 (right). Its real 
(dashed line) and imaginary (dotted line) parts are compared to the harmonic oscillator (solid line). 



whose general solution, for any e G C, is given by 

1-e 1 



u(x) = exp(— x /2) 



4 '2' 



+ (A + u;)xiFi ( — - -, ^;x 2 



(13) 



With this formalism, the first-order SUSY partner potential V of the harmonic oscillator is 
given by 

V{x) =x 2 -2[lnu(x)]". 

The previously known results for the real case [27] are obtained by taking e E K, k = and 
expressing A as 



A = 2v 



r(¥) 
r(¥)' 



On the other hand, for e G C the transformation function u(x) is complex and so is V{x). 

In Fig. 2 we present some examples of complex SUSY partner potentials of the harmonic 
oscillator generated for e G C and compare them with the initial potential. Let us note that 
these new potentials have the same real spectra as the harmonic oscillator, except that they have 
one extra energy level, located at the complex value e. This kind of spectrum is represented in 
the diagram of Fig. 3. Note that this can be interpreted as the superposition of the two ladders 
shown in equation (7). 



5 Complex solutions to Painleve IV equation 

In this Section we will see that for each first-order SUSY partner of the harmonic oscillator we 
will have three different complex solutions of Pjy. 

First of all, let us note that the ladder operators associated with H are given by 

L ± = A + a ± A~, 

or, explicitly, 



l+ = (-d + p)(-d + x)(d + (3) = (d - /3)(d - x){d + p), 
l- = (-d + p){d + x)(d + p) = (-d + p)(-d - x)(-d - p). 



(14) 
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Figure 3. The complex energy plane which contains the eigenvalues of the SUSY generated Hamil- 
tonian H . Along the real line it is drawn the usual ladder formed by the energy levels E n = 2n + 1, 
n = 0, 1, ... , and out of the real line there is a new level of the complex value e. 



This system has third-order differential ladder operators, therefore it is ruled by a second-order 
PHA and so we can apply our analysis of Section 3 by identifying the ladder operators given 
in equations (14) with those of equations (8). This identification leads to 

fi = -P, f2 = x, f 3 = /3. (15) 

Recall that the function g = — x fulfills the Painleve IV equation, then one solution to F\y is 

g = (3 - x. 

Let i/jg i , i = 1,2,3, be the states annihilated by L~ , where Ei represents the factorization 
energy for the corresponding extremal state. In particular, let tp£ 3 be annihilated by LJ, and 
from equation (8) we can see that it is also annihilated by L~ . By solving L^ip£ 3 = we get 



tp£ 3 oc exp 



h{y)dy 



from which it can be shown that the corresponding solution to Pry reads 
93(x,e) = -x - ]n[if) g3 (x)]' . 

For our complex first-order SUSY partner potential V(x), generated by using the complex 
seed solution u(x) of equation (13), the three extremal states (up to a numerical factor) and 
their corresponding energies consistent with equations (15) are given by: 



ip £l = A + a + u, £x = e + 2, Vf 2 = A + exp (-x 2 /2) , £ 2 = 1, 



(16) 



Nevertheless, the labelling given in equations (16) is not essential, so that by making cyclic 
permutations, associated with the three extremal states of equations (16) we get three solutions 
to Piv: 



9i(x,e) = -x - ln[ip Si (x)]' , 
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-4 -2 2 4 



Figure 4. Complex solutions to Pry- The solid line corresponds to the real part and the dashed to the 
complex one. gi(x) for e = —1 + il0~ 2 and A = k = 1; 32(2;) for e = 4 + «2 _1 , A = ac = 1; and gs(x) is 
given for e = 1 + i and A = 3, k = 1. 

with i = 1,2, 3. The corresponding parameters a, 6 of Pry are given by: 

1 1 2 

ai = --(e + 5), &i = --(e-l) , a 2 = e - 1, 6 2 = -2, 

a 3 = ^(l-e), 6 3 = -^( e + l) 2 , (17) 

where we have added the subscript corresponding to the extremal state used. In Fig. 4 we have 
presented one example for each of the three families of solutions. 

From equations (17) we can see that Oj is linear in e for all three cases. So, instead of studying 
the parametric relationship of a, b in terms of e let us analyze b{ = bi(di), namely, 

fei = -2(a!+3) 2 , 62 = -2, 63 = -2(a 3 - l) 2 . (18) 

Then, we can choose ai G C but 6, will be fixed by its corresponding relationship with a\. In 
Fig. 5 we show the domain for 63 from equations (18). There are similar plots for b\ but for 62 
remember that 62 = — 2 Va2 G C. 

6 Conclusions 

In this work we have applied the first-order SUSY QM to the harmonic oscillator potential to 
obtain systems with third-order ladder operators which are ruled by second-order PHA. From 
this system we have obtained complex solutions g\, g2, and 53 of the Painleve IV equation. 

It is worth to notice that complex solutions have not been generated in this context, never- 
theless, the work of Bassom et al. [8, Section 3.3] is general and can be applied to a Pry with 
complex parameters. An extension of the parameter space to the complex plane is helpful at 
least when P\y appears in the reduction of other non-linear differential equations. Moreover, it 
can be used in the future to study complex extensions of those physical systems described by P\y. 
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Figure 5. Parameter space where we show Re(63) (left) and Im(&3) (right) in function of Re(<23) and 
Im(d3). The plots for b\{a\) and 62(02) are similar. 

We first discussed the complex first-order SUSY transformation for general systems. We 
defined then the PHA and established explicitly the relationship of the second-order case with 
the Painleve IV equation. We also studied the complex SUSY partner potentials of the harmonic 
oscillator with an analysis of the energy spectra of the new potentials. Then, we obtained 
complex solutions of the Painleve IV equation and we analized the complex parameter space 
where the solutions can appear. We concluded that there are three different families of complex 
solutions for each complex first-order SUSY transformation. 

In the future we would like to address the higher-order SUSY transformation in order to 
extend the space of complex solutions of the Painleve IV equation and to study in more detail 
the properties of these solutions. 
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